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Abstract 



Q, 

, We present a pilot narrow-band survey of Ha emitters at z = 2.2 in the Great Observatories Origins 

■ Deep Survey North (GOODS-N) field vifith MOIRCS instrument on the Subaru telescope. The survey 
reached a Scr limiting magnitude of 23.6 (NB209) which corresponds to a 3ct limiting line flux of 2.5 x 

. 10~^^ erg s^^ cm^^ over a 56 arcmnin^ contiguous area (excluding a shallower area). From this survey, 

' we have identified 11 Ha emitters and one AGN at z = 2.2 on the basis of narrow-band excesses and 

, photometric redshifts. We obtained spectra for seven new objects among them, including one AGN, and 



> 



an emission line above Scr is detected from all of them. We have estimated star formation rates (SFR) and 
Q^l stellar masses (Mstai) for individual galaxies. The average SFR and Mgtar is 27.8 MQyr"-"^ and 4.0x 10^°Mq, 

I ' respectivly. Their specific star formation rates are inversely correlated with their stellar masses. Fitting to 

Q I a Schechter function yields the Ha luminosity function with logL = 42.82, logcj) = —2.78 and a = —1.37. 

The average star formation rate density in the survey volume is estimated to be 0.31 Moyr~^Mpc^'^ 
' according to the Kennicutt relation between Ha luminosity and star formation rate. We compare our Ha 

emitters at z = 2.2 in GOODS-N with narrow-band line emitters in other field and clusters to see their time 
evolution and environmental dependence. We find that the star formation activity is reduced rapidly from 
z=2.5 to 2=0.8 in the cluster environment, while it is only moderately changed in the field environment. 
' This result suggests that the timescale of galaxy formation is different among different environments, and 

\^ I the star forming activities in high density regions eventually overtake those in lower density regions as a 

00 ■ consequence of "galaxy formation bias" at high redshifts. 

' Key words: galaxies: evolution - galaxies: formation - cosmology: early universe 

1. Introduction (Kodama et al., 1998). The blue star- forming galaxies, 

however, are much more difficult to be identified at z ~ 2 
Since recent observations in optical and near-infrared using a color-selection or a photometric redshift technique 
_ ^ wavelength indicate that the volume-averaged star forma- because their spectra are relatively flat and featureless in 
, tion rate (SFR) increases from z = to z ^ 1 and plateaus the optical-NIR regime. Moreover, in a narrow-band sur- 
' at z ^ 2 (Hopkins, 2004; Hopkins & Beacom, 2006), it is vey which captures emission lines from galaxies directly, 
5^ , likely that a large fraction of stars in galaxies at present- we can sample star forming galaxies completely above a 
day formed at z > 1. The AGN activity and the red- certain limiting flux and an equivalent width limit unless 
shift distribution of submm galaxies (SMG) also peak at the line is attenuated by dust or stellar absorption. We 
this epoch (Chapman et al., 2005). Therefore the redshift are not biased by colors of galaxies, either. For these rea- 
rangc of z=2-3 is the epoch when galaxies have the most sons, the narrow-band emitter survey is a very efficient 
intensive evolution. It is absolutely imperative to build and effective method for initially making a sample of star 
a statistical sample of star forming galaxies at z=2-3 in forming galaxies and obtaining their photometric proper- 
order to understand the cosmic star formation history and ties at a particular redshift. We can then study them more 
the early evolution of galaxies. A line-emitter survey with in detail by spectroscopic follow-up in the near-infrared 
a narrow-band filter can provide a large sample of active and also by radio observations targeting their molecular 
galaxies efficiently in a limited range of redshift in contrast lines. 

to color-color selections such as UGTZ (Stcidel et al., 2004) One of the most important quantities characterizing a 

or BzK (Daddi et al., 2004). Since red galaxies represent galaxy is SFR. There are many different SFR indicators, 

characteristic colors on the color-magnitude diagram as such as UV continuum, nebular emission lines (Ha and 

referred to as the "red sequence", we can relatively eas- [OH]) and mid- infrared, but the SFRs estimated by vari- 

ily select such red galaxies located at a specific redshift ous measurements do not always provide consistent results 
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due to selection biases and different amounts of dust ex- 
tinction effects, and so on. For UV continuum radiated 
by hot OB type stars, we must correct for absorption by 
the surrounding dust. This effect is often corrected for by 
using the UV slope but this process may lead to a large 
uncertainty in the estimated SFRs. On the other hand, 
mid-infrared emission is not affected by dust "extinction" 
since it is dust "emission" . However, we can observe only 
the galaxies that have relatively high SFRs so as to be 
detected by mid-infrared observations which are less sen- 
sitive to optical or near-infrared observations. Also, to 
derive SFRs from the mid-infrared luminosity one has to 
assume the dust temperature which is somewhat uncer- 
tain. The Ha line, a hydrogen's Balmer series line emit- 
ted from ionized gas (i.e., HH region) around hot young 
stars, is one of the best SFR indicators. It has many 
great advantages; being less affected by dust extinction, 
providing a survey with high sensitivity, and having been 
well calibrated in the local Universe. For the rest-frame 
UV-selected galaxies, Erb et al. (2006a) obtained 114 Ha 
spectra of star forming galaxies at z ~ 2 and indicate that 
Ha emission is attenuated by a typical factor of 1.7, which 
is about 1/3 compared to the UV attenuation. 

The narrow-band survey of Ha line enables us to make 
a relatively-unbiased large sample of star forming galaxies 
at the same redshift as well as to obtain accurate estimates 
of SFR, unlike the rest-frame UV surveys which tend to 
miss a significant fraction of star-forming activities due 
to dust obscuration. However, the Ha line is red-shifted 
into the near-infrared regime at z > 0.5. Because the near- 
infrared observation is severely affected by OH sky emis- 
sion lines, only a fraction of redshift range can be surveyed 
with high sensitivities. At z ^2.2, Ha line falls just in be- 
tween the OH line forest (A ^ 2.1^m). There are some 
pioneering previous works that searched for Ha emitters 
at z^2 with narrow-band imaging (Thompson et al., 1994; 
Bunker et al., 1995; Teplitz et al., 1998) or with slitless 
spectroscopy (Yan et al., 1999; Atek et al., 2010). In the 
past decade, some large Ha surveys have been carried 
out at z ^ 2.2. Moorwood et al. (2000) searched for Ha 
emitting galaxies with ESO NTT telescope. The survey 
reached a limiting line flux of 5xl0~^^ erg s~^ cm~^ and 
covered a 100 arcmin^ area. The observation with the 
narrow-band filter yielded 10 candidates and 6 of them 
have been confirmed spectroscopically later. Then, Geach 
et al. (2008) conducted a narrow-band survey with a line 
fiux limit of ^ 10~^^ erg s""'^ cm~^ over a large area of 
0.6 deg^ with a 3.8m telescope (UKIRT) and identified 
55 Ha emitters. Although their survey provided a sample 
large enough to evaluate the Ha luminosity function (LF) 
at the bright end, they pointed out that the faint end of 
Ha LF is important in order to estimate the global star 
formation rate, and hence a deeper survey is required for 
this purpose. To date, the deepest survey of Ha emitters 
at this high redshift is that of Hayes et al. (2010). They 
carried out a narrow-band survey with HAWK-I on ESO- 
VLT, and reached to a 3cr flux limit of 4.1 xlO""'^^ erg 
s~^ cm~^ over 56 arcmin^ area. Their survey suggested 
that the faint-end slope of the Ha LF is steeper than the 



value obtained in the local universe. However, large scale 
structures at z ^ 2 may stretch the intrinsic slope of the 
LF and it is too early to conclude on the faint end of Ha 
LF by only the data in a single fleld. The result is very 
susceptible to the effect of cosmic variance. 

Although the recent development of wide-field instru- 
ments in the near-infrared regime has enabled us to study 
distant galaxies at z ~ 2 systematically, there is an on- 
going debate as to how the ancient galaxies evolve into 
the ones at the present-day universe. This is attributed to 
the fact that the galaxy evolution is connected not only 
to time, but also to environment and mass (e.g. Elbaz 
et al., 2007; Tanaka et al., 2010). It is thought that there 
are two types of causes for the environmental dependence, 
inherent origin and acquired effects. In the former case, 
initial density fluctuation is originally a bit large in high 
density regions such as cluster cores, and galaxy forma- 
tion take place earlier there compared to lower density 
regions. In the later case, some external factors, such as 
gas stripping by ram pressure and mergers between galax- 
ies, influence galaxy properties during the course of their 
build-up. The spatial distribution of star forming galaxies 
helps us to tell whether the inherent origin or the acquired 
effects dominant the environmental dependence. In galaxy 
cluster RXJ1716 at z = 0.8, Ha emitters or mid-infrared 
sources, which are both tracers of active star formation, 
are preferentially found in the surrounding region of the 
cluster core (Koyama et al., 2008, 2010). In contrast, star- 
forming galaxies in XCS2215 at z = 1.5 are located in the 
cluster core as often as in the surrounding environments 
(Hayashi et al., 2010; Tran et al., 2010). It seems that 
star formation activity is initially enhanced in massive 
galaxies in high density regions and is propagated to lower 
mass galaxies and to lower density regions with time. If 
the relative velocity between galaxies is too large, these 
galaxies would pass through each other without merging 
together. Therefore, at z = 0.8, it is thought that mergers 
in mid-density regions would be more important than the 
inherent effect. However, at a higher redshift z = 1.5, it 
seems that the inherent effect would surpass the acquired 
effects. To further look into and confirm this interesting 
result based on the two clusters, and to extend it further 
back in time, it is required to construct a larger sample 
of star forming galaxies over the critical era of z =1-3, 
that covers various environments from galaxy cluster core 
to blank field region in order to compare the properties 
between different environments. 

In this paper, we present a survey of Ha emitters at 
z = 2.2 in GOODS-N field in order to study star formation 
activities in the general field or low density regions. This 
paper is structured as follows. In § 2, we describe our sur- 
vey design, observations, and the data. In § 3, our target 
selection of Ha emitters and the spectroscopic follow-up 
observations are described. We show our results in § 4 
and § 5, and discuss the star formation activities of galax- 
ies in various environments and redshifts in § 6. Finally, 
we summarize our study in § 7. Wc assume the cosmo- 
logical parameters of Ho=70 km s~^ Mpc~^, JIm = 0.3, 
and Oa = 0.7, and adopt AB magnitudes throughout this 
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Fig. 1. The transmission curves of our triple narrow-band 
filter system (black) and broad-band filters (red) . Lower panel 
shows the wavelength of OH sky emission lines (Rousselot 
et al., 2000). Blue solid line shows the template spectrum of 
a typical star-forming galaxy at z = 2.2 (Kinney et al., 1996). 



paper. 
2. Data 

2.1. Triple narrow-band filter system for z ~ 2.2 

Taking the advantage that Ha Hne at z = 2.19 is free 
from any strong OH sky emission hne, we have made 
three sets of custom-made narrow-band fihers, namely, 
NB119 (Ac = 1.19^m, FWHM=0.014^m), NB155 (Ac = 
1.55^m, FWHM=0.017/im) and NB209 (Ac = 2.09^m, 
FWHM=0.027^m). Figure 1 shows the filter response 
functions of these narrow-band filters together with those 
of broad-band fihers J, H and K^. NB119, NB155 and 
NB209 filters capture [OH], H/3 and Ha emission lines, 
respectively, from the same galaxies at z = 2.19 ±0.02. 
If emission lines are detected in more than a single filter 
from a galaxy, we can confirm its redshift without spec- 
troscopy. Even if only a single line is detected, all we need 
to do is to discriminate among [OH]/[OHI]/Ha by using 
photometric redshifts. The combination of NB209/NB155 
filters tries to measure Balmer decrement (Ha/H/3) in or- 
der to accurately correct for dust extinction for the bright- 
est emitters. In this paper, however, we will concentrate 
on the studies using the NB209 filter only, since we have 
not obtained deep enough data with NB155 yet. 

2.2. Imaging data 

The imaging survey has been carried out with MOIRCS 
(Suzuki et al., 2008) on the Subaru telescope (lye et al., 
2004) with the NB209 fiher in GOODS-N field. MOIRCS 
is equipped with two Hawaii-2 detectors (2048x2048). 
Unfortunately, however, because one of the two MOIRCS 
chips was unavailable at the time of the observation run 
due to a trouble, one pointing could cover only a half of 
the original FoV that is 4' x 3.5'. We spent five point- 
ings in total in April 2008 and the total area of the sur- 
vey amounts to 70 arcmin^, corresponding to about 9300 
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Fig. 2. The K - NB209 color-magnitude diagram to select 
Ha emitter candidates at z = 2.2. Grey dots show all the 
objects detected in both Kb and NB209 images and black 
dots show the NB209 emitters, which satisfy the criteria 
in section 3.1. Ho emitters and [OIH] emitters are shown 
by blue and red squares, respectively. The solid curve 
and the dashed line show 3cr photometric error and nar- 
row-band excess of 0.4 magnitude in Ks - NB209, respectively. 



Mpc'^ in co-moving volume for the narrow-band survey 
at z = 2.19 ±0.02. Table 1 shows the observed fields, 
exposure times, point spread functions (PSFs) and lim- 
iting magnitudes of NB209 and ii's-band. Data reduction 
was conducted by using MCSRED (I. Tanaka, et al., in 
preparation). Selfsky flat image was used to correct for 
variabilities in sensitivity from pixel to pixel. Before com- 
bining all the images, PSFs were matched to 0.8" in GT4 
and to 0.7" in other fields by gaussian smoothing. We 
calibrated the flux by observing a spectroscopic standard 
star, G191-B2B, and measured the photometric zero-point 
of our NB209 image (23.9). Limiting magnitude was esti- 
mated by using "limitmag" , which is a task of SDFRED 
(Yagi et al., 2002; Ouchi et al., 2004). Note that the 
narrow-band image in GT4 is shallower than the other 
flelds. To select line emitters, and to measure their line 
fluxes, the continuum levels should be measured as well. 
As shown in Fig. 1, the NB209 filter is in the middle of the 
Ks filter and Kg magnitudes can be used to measure the 
continuum levels. For Kg band data, we used the data- 
set of the MOIRCS Deep Survey (MODS; Kajisawa et al. 
2009), which covers 70% of the GOODS-N field. PSF of 
the Ks band data was matched to the NB209 data. Also, 
we used the MODS photometric catalog (Kajisawa et al., 
in prep), which contains optical (U BVi' z'), near-infrared 
(JHKs), mid-infrared (IRAC 3.6, 4.5, 5.8, S/im and MIPS 
24/xm) and X-ray data. The photometric redshifts and the 
stellar masses, based on GALAXEV (Bruzual & Chariot, 
2003), are also obtained from this catalog. 
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Table 1. Summary of the observations. 
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Table 2. Summary of our 12 H« emitters. 



ID 




a 

[J2000] 


8 

[J2000] 


A's 










/"line 

[10- 


17 


erg s ^ cm ^ ] 


SFR 

[Moyr- 




■^spcc 


comment 


MODS- 


91 


189.0141 


62.1760 


23 


07 


± 





16 


2.28 


± 


1.82 


12.2 ± 


9.7 








MODS- 


487 


188.9930 


62.1988 


23 


19 


± 





17 


3.57 


± 


2.06 


18.8 ± 


10.8 


2 


176 




MODS- 


3653 


189.0906 


62.2480 


22 


65 


± 





13 


3.09 


± 


2.25 


16.7 ± 


12.2 


2 


200 




MODS- 


6422 


189.2159 


62.2513 


22 


34 


± 





11 












2 


193 


AGN 


MODS- 


7651 


189.2582 


62.2639 


23 


37 


± 





18 


5.93 


± 


2.22 


31.6 ± 


11.8 


2 


188 




MODS- 


7773 


189.2364 


62.2788 


22 


53 


± 





13 


5.43 


± 


2.39 


29.0 ± 


12.8 


2 


190 




MODS- 


7889 


189.2204 


62.2907 


23 


53 


± 





20 


3.87 


± 


1.91 


20.6 ± 


10.2 


2 


187 




MODS- 


7899 


189.2235 


62.2901 


21 


94 


± 





09 


11.2 


± 


3.32 


59.8 ± 


17.7 


2 


186 




MODS- 


8313 


189.2152 


62.3071 


24 


00 


± 





25 


3.98 


± 


1.71 


21.2 ± 


9.1 








MODS- 


8540 


189.2709 


62.2908 


22 


63 


± 





13 


8.54 


± 


2.83 


45.6 ± 


15.1 








MODS- 


9441 


189.3259 


62.2815 


22 


43 


± 





12 


4.26 


± 


2.46 


22.9 ± 


13.2 


2 


196 




MODS- 


10522 


189.3564 


62.3194 


22 


95 


± 





15 


5.17 


± 


2.16 


27.7 ± 


11.6 


2 


192 





use A's image for a detection. It is therefore unlikely that 
we miss many high equivalent width emitters with low 
continua. The extraction criteria were at least 25 pixels 
with fluxes above la level, where Icr is the sky noise of the 
image. (2) Marginal objects, whose fluxes were less than 
5cr in NB209 or less than 3(T in A's band, or those located 
at the edge of the image, were rejected as fakes. (3) In 
order to avoid statistically scattered objects around zero 
in /\s-NB209 color due to photometric errors, we selected 
only the objects whose excesses were larger than 3cr or 
0.4 magnitude, corresponds to about 40 A in the equiva- 
lent width in the rest frame (Fig. 2). Finally, 21 emitter 
candidates were selected. 

It is expected that the NB209 filter detects not only 
Ha emitters at z = 2.2 but other line emitters at different 
redshifts. For example, there may be contamination of 
[OIII] emitters at z = 3.2 and [Oil] emitters at 2; = 4.6. In 
order to separate Ha emitters from the others, we used 
the spectroscopic redshifts and the photometric redshifts, 
taken from the MODS catalog. Four out of 21 candidates 
had spectroscopic redshifts, i.e. 2=0.840, 2.186, 2.200 and 
3.187. Figure 3 shows the distribution of the photometric 
redshifts for the NB209 emitters. This clearly exhibits a 
bimodal distribution at z = 2.2 and z = 3.3. The criterion 
of 1.8 < Zphot < 2.4 was adopted to select Ha emitters, and 
we identified 12 Ha emitters. All of them also satisfy the 
BzK criterion (Daddi et al., 2004). One of them is likely 
to be a AGN because a X-ray emission is detected from 
this object (Alexander et al., 2003). The coordinates and 
the properties of these final emitters are listed in Table 2. 
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Fig. 3. The distribution of photometric redshifts for 
the NB209 emitters. The vertical dashed lines show 
our criteria to select Ha emitters (1.8 < Zpijot < 2.4). 



3. Target selection 

3.1. Narrow-band emitters 

Emission line objects were selected on the basis of the 
excesses of NB209 fluxes over Ks fluxes with the following 
criteria and procedures: (1) Source detection was per- 
formed on the Ks band image using SExtractor (Berlin & 
Arnouts, 1996). Usually, a narrow-band image is used for 
a detection. However, since our Kg band image is much 
deeper than the NB209 image by 1.5-2 magnitudes, we 
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3.2. Spectroscopic confirmation 



servation in the J-band. 



For the seven new emitters out of 12, including the 
AGN candidate, we obtained near-infrared spectra with 
MOIRCS on multi-object spectroscopy (MOS) mode 
(Tokoku et al., 2006). Our aims are to confirm that 
the narrow-band emitters are the Ha emitters located at 
z = 2.2, and to accurately measure their redshifts. The 
medium-resolution grizm, R1300-I-K, was used to possibly 
discriminate between Ha and [NH] lines. The slit width 
was 0.7" and its length was 10-12". The on-source expo- 
sure time was 320 minutes. Data reduction was conducted 
with the following procedures. First of all, a Sky frame 
was created by combining all the object frames. We define 
sky noise as the poisson noise of this frame at each wave- 
length. The Sky frame and the object frames were then 
divided by the dome fiat for fiat fielding. Cosmic-rays and 
bad pixels were removed. Distortion was corrected by us- 
ing the task "mcsgeocorr" of MCSRED. The wavelength 
calibration was performed by using OH lines in the Sky 
frame. For object frames, the Sky pattern was roughly 
removed by subtracting the adjacent dithered image. The 
residual background emission was subtracted by using the 
task "background" of IRAF. To calibrate the telluric ab- 
sorption and the difference in sensitivity at each wave- 
length, the object frames were divided by the spectrum of 
a spectroscopic standard star (AO- type), and multiplied 
by the model spectrum of a AO- type star (Pickles, 1998). 
Flux is calibrated by using a star in the 2-MASS catalog 
(Skrutskie et al, 2006). 

Figure 4 shows their spectra in the rest frame. An emis- 
sion line, whose flux is above Scr, is detected from all of 
them. For MODS-7773, we can also see a weak [NH]A6584 
line, hence these objects must be Ha emitters for sure. For 
the others, [NH] lines do not appear. However, if a NB209 
emitter is a [OHI] emitter at z—3.2, the doublet lines of 
[OHI]AA4959,5007 (in the rest frame) or Aobs=2.083, 2.103 
(in the observed frame) should have been resolved, given 
the resolution of our observations. But we do not detect 
any. Moreover, since the ratio of the [OHI] doublet is 
1/3 (Arnaboldiet al., 2003), we can reject the possibility 
of [OHI] emitters at z = 3.2 for MODS-7651, 9441 and 
10522, whose lines are detected at more than 5cr, and the 
weaker line should have been visible. If any of them is a 
[OH] line at z = 4.6, SFR based on the [OH] fiux, would 
be estimated to be about 200 Moyr^^. Despite of such 
very high SFRs, however, these objects are not detected 
by MIPS 2Afim {S24 < 80/xJy). Therefore, these objects 
are also likely to be Ha emitters at z = 2.2. MODS-6422 
is thought to be an AGN due to the X-ray detection. The 
Ha profile looks marginally extended and the FWHM of 
the velocity width is about 850km s~^, which is more than 
twice wider than those of other Ha emitters. It is thus con- 
firmed that our criteria, satisfying the NB209 excess and 
1.8 < Zphot < 2.4, is effective and robust in removing con- 
tamination and selecting genuine Ha emitters at z ~ 2.2. 
In order to make more accurate confirmation, however, 
one would need to detect [OH] lines, which redshift to 
1.19/im at z = 2.2, by another spectroscopic follow-up ob- 



4. Star formation rate and stellar mass 

We have constructed a sample of 11 Ha emitting star- 
forming galaxies at z = 2.2 in the GOODS-N field. We 
now estimate their SFRs and stellar masses. From the 
flux densities in NB209 and Kg band, i.e. /nb and /bb, 
respectively. Ha line fluxes can be calculated as follows, 

AnbAbb 



-^linc 



iBB 



iNB 



-(/nb - /bb), 



(1) 



where A indicates a wavelength width (FWHM) of 
a filter and / denotes the fiux density. We use 
ANB=0.027[/xm] and ABB=0.31[/xm]. Wc note that the 
[NH] line also contributes to a NB209 flux. We adopted 
L([NII])/L(Ha-|-[NII])=0.25 for afl the objects (Sobral 
et al., 2009). Using the standard calibration of Kennicutt 
(1998a); SFR[Moyr-i]=7.9 x lO^''^ L(Ha)[erg s'^], we 
converted the Ha luminosity to a SFR, assuming the 
Salpeter IMF (Salpeter, 1955). To correct for dust extinc- 
tion, we derived E{B — F)stciiar from the SED fitting. In 
the local universe, it is known that emission lines from the 
ionized gas are more attenuated than stellar continua such 
as E{B- F)steiiar = O.UE{B - F)gas (Calzctti et al. , 2000) . 
On the other hand, at z ~ 2, Erb et al. (2006) used the 
same extinction values for both stellar continua and the 
nebular emission lines, i.e. E{B — F)steiiar = E{B — T^)gas- 
This is because the above formula by Calzetti et al. (2000) 
significantly overestimates the Ha-based SFRs with re- 
spect to the UV-based SFRs. The average strength of 
dust extinction for our Ha emitters is estimated to be 
A(Ha)=0.7-1.6 if the Calzctti formula is used. In this pa- 
per, we adopt the constant value of A(Ha)=l, given the 
uncertainty. As a future work, however, the amount of 
reddening should be ideally estimated from the Balmer 
lines ratio, such as Ha/H/3. For emitters with spectro- 
scopic redshifts, we also correct for the filter transmission 
depending on the location of the Ha line on the response 
function of the narrow-band filter. For those without 
spectroscopic redshifts, we assumed the average redshift 
of z = 2.190. The stellar masses were obtained from the 
MODS catalog (Kajisawa et al., 2009) which were derived 
from the SED fitting. The average SFR and Afgtar of our 
Ha emitters are 27.8 Mq yr~^ and 4.0x lO-^^Af©, respec- 
tively. 

The specific star formation rate (SSFR), defined as the 
SFR per unit stellar mass, shows how much the current 
star formation activity contributes to increase the stellar 
mass of a galaxy. For UV and near- infrared selected galax- 
ies, it is reported that there is a good correlation between 
SSFR and stellar mass over a wide range in the evolution- 
ary stage from young galaxies which begin to form stars, 
to passively evolving galaxies (Erb ct al., 2006a; Reddy 
et al., 2006). SSFR is high in a less massive galaxy and 
it decreases with increasing stellar mass. Figure 5 shows 
SSFR versus stellar mass for our Ha emitters. Open cir- 
cles show the UV-selected galaxies of Erb et al. (2006a, b). 
Because they adopted the Chabrier IMF (Chabrier, 2003), 
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we convert their data points to those for the Salpeter IMF 
by multiplying 1.8 to the SFRs. For our Ha emitters, the 
SSFR negatively correlates roughly with the stellar mass 
as seen in the samples of Erb et al. (2006a, b). Our sam- 
ples have systematically lower SSFRs than those of the 
UV-selccted samples. Two selection effects would be con- 
tributing to this offset. One is that our survey is deeper 
and can detect galaxies with much weaker SFRs down to 
12 Mgyr-^ In other words, it may merely reflect the dif- 
ference in the flux limits. The other is the large difference 
in the survey volumes. Because the survey area of Erb et 
al. (2006a,b) is very wide, 0.48 deg^ (Steidel et al., 2004), 
corresponding to 6.1 x 10^ Mpc^ in comoving volume at 
z = 1.4 — 2.5, their survey can efficiently detect relatively 
rare objects with strong SFRs (e.g., SFR - 100 Mgyr-^). 
Based on the luminosity function derived in § 5, we can 
expect to detect only 2 high SFR galaxies (> 100 MQyr~^) 
m our survey volume (9300 Mpc^). Therefore, the non- 
detection can be simply due to a cosmic variance effect. 
The inverse of SSFR indicates the star formation time- 
scale to form all the stars in a galaxy. The average time- 
scale of the Ha emitters is (r) = 1.4 Gyr, which is smaller 
than the age of the universe at z—2.2 (3 Gyrs). Therefore 
star formation history of our Ha emitters does not neces- 
sarily require a star-burst and can be reproduced even by 
a constant star formation model. 




10'" 10' 

Mstar [Mq ] 



Fig. 5. The specific star formation rates (SSFR) versus stel- 
lar masses. The red diamonds and the open circles in- 
dicate our Ho emitters and the UV-selected galaxies (Erb 
et al., 2006a, b), respectively. We assumed the dust extinc- 
tion of A(Ha)=l [mag] and Salpeter IMF for both samples. 
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5. Luminosity function and cosmic star forma- 
tion rate density 

In order to estimate tlie volume-averaged star formation 
rate density, it is necessary to extrapolate the luminos- 
ity function to the unobserved faint star-forming galax- 
ies. The Ha LF is well represented by a Schechter func- 
tion (Schechter, 1976). In the local universe, the faint- 
end slope is measured to a = —1.3 (Gallego et al., 1995). 
However, recent studies indicate steeper faint-end slopes 
(a ^ -1.6) in the distant universe (Hayes et al., 2010). 
Our deep survey can constrain the faint-end slope of Ha 
luminosity function, independently. The volume density 
of galaxies is computed as 

* = ir^y^, (2) 

AlogL^14iax ^ ' 

within a luminosity bin of A log L = 0.4. Vinax is the max- 
imum co-moving volume probed by the survey. Because 
the filter response function does not have a perfect rect- 
angle shape, Vinax is larger for more luminous galaxies. 
We have calculated the luminosity function excluding the 
GT4 field because the depth of the narrow-band images 
in this filed is significantly shallower than the other fields. 

Figure 6 shows the Ha luminosity function at z ~ 2.2 
obtained by this survey. Ha luminosity has not been cor- 
rected for dust extinction. For the bright end. Ha emit- 
ters of Geach et al. (2008) are used. Fitting a Schechter 
function to the data, without any correction for dust ex- 
tinction, gives logL = 42.82 ±0.65, log(?!)= -2.78 ±1.08 and 
a = —1.37 ±0.44. Our result does not support the steep 
slope at the faint end. This indicates that the cosmic vari- 
ance has a considerably large impact on the result. Figure 
7 shows the spatial distribution of Ha emitters identified 



by this survey. Clearly, they are not homogeneously dis- 
tributed across the field. This is not explained by the 
difference in the depth of the imaging data among differ- 
ent fields (Table 1). Most of our samples are concentrated 
in the GT3 field and there are only a few emitters in the 
other fields. We have calculated a LF in GT3 field only to 
investigate the influence of the cosmic variance. Fitting 
a Schechter function yields logL = 43.36, log0 — —3.70 
and a = —2.08, which now supports the steep slope at the 
faint end. This suggests that it is difficult to derive a rep- 
resentative number density of Ha emitters at z = 2.2 by a 
small survey (<100 arcmin^). In contrast, the survey area 
0.6 deg^) of Geach et al. (2008) is wide enough to av- 
erage the bright end of LF over the large-scale structures 
at z = 2.2. It is thus also required to have an additional 
deep survey (Lua ^ 10*^[erg s~^]) over a much wider area 
to negate the effect of cosmic variance. 

From the above derived Ha LF, we can estimate the 
cosmic star formation rate density (SFRD). The Ha LF is 
integrated over the range of 10'^^ <LHQ[erg s ~^]< 10^^ as 
in Hopkins (2004). The resultant volume-averaged SFRD 
is PsFR = 0.31 MQyr~^Mpc~'^. If it is integrated only 
over the range of 10^^ <LHct[erg s ~^]< 10^^, the resultant 
SFRD is 0.21 Mgyr-iMpc-^. If the faint end slope of the 
LF is steeper, the proportion of the extrapolation at the 
faint end would increase. The derived SFRD would be 
nearly doubled in the case of a = —1.6. The evolution of 
SFRD has been studied with various indicators. We com- 
pare our result with those of previous studies with Ha 
narrow-band imaging, and show the evolution of SFRD in 
Figure 8 (Perez-Gonzalez et al., 2003; Dale et al., 2008; 
Morioka et al, 2008; Ly et al., 2007; Shioya et al., 2008; 
Fujita et al, 2003; Hippclcin et al., 2003; Vfllar et al, 
2008; Sobral et al., 2009; Moorwood et al., 2000; Geach 
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et al., 2008; Hayes et al., 2010). The SFRD increases 
up to z ^ 1 but there is little change between z = 1 and 
z = 2. We note that this trend is also consistent with 
the result based on the spectroscopic sample at 2:=0.8, 
that is 0.1 Mgyr-^Mpc-^ (Dohcrty et al., 2004, 2006), 
although the direct comparison is difficult due to differ- 
ences in sample selection and completeness. According to 
other star formation indicators, the SFRD does not in- 
crease but rather decreases at z> 3 (Hopkins &: Beacom, 
2006). Because Ha shifts to A > 2.6/im at z > 3, it is diffi- 
cult to detect Ha emission lines from distant galaxies with 
the ground-based telescope. The space infrared telescopes 
such as JWST and SPICA (Nakagawa, 2008) will be very 
powerful for tracking the precise star formation histories 
further back in time at z > 3. 

6. Star formation history as functions of redshift 
and environment 

As mentioned in § 1, the evolution of galaxies depends 
strongly on the environment. Because field galaxies are 
antithetical to cluster galaxies in terms of environment, 
we can directly examine the environmental dependence 
of star formation activity by comparing the properties of 
emitters among different environments. In addition to the 
present general field survey, we have estimated SFRs and 
stellar masses for a large sample of star forming galaxies 
both in clusters and in the other general field at various 
redshifts (Field: [OH] emitters at z = 1.2 in COSMOS, 
Ideue et al. 2009; Clusters: Ha emitters at z — 0.8 in 
RXJ1716, Koyama et al. 2010; [OH] emitters at z = 1.5 
in XCS2215, Hayashi et al. 2010). All the samples are se- 
lected on the basis of narrow-band nebular emission line 
surveys (Ha or [OH]). We have limited the galaxy sam- 



ple to those having SFRs above 12 MQyr^^ and stellar 
masses above lO^*' M0 to make fair comparisons. We 
have estimated the stellar masses of the emitters with us- 
ing the theoretical mass-to-luminosity ratios predicted by 
the population synthesis model (Kodama et al., 1999) af- 
ter it is scaled to the Salpeter IMF for consistency. At 
z = 2.2. 



log(Af,tar/10") = -0.4(/^f* - 

AlogM = 0.229[(z - A's) - 3.45]. 



(3) 
(4) 



z and Ks indicate z-band and ii's-band magnitudes, re- 
spectively. K^^ denotes the A'g-band magnitude for a 
galaxy that has a stellar mass of 10^^ Mq, and it is esti- 
mated to be 22.42 at z = 2.2. The upper equation exhibits 
the mass-magnitude relation in the case of a passively 
evolving galaxy. The lower equation gives an amount of 
correction for stellar mass depending on the z — K color, 
with respect to the mass given by the upper equation. The 
mass-to-light ratios of galaxies are different from redshift 
to redshift and among different colors, and are estimated 
by the model. The Ha/ [OH] ratio is assumed to be log 
(Ha/[OII])obs=0.2 which is the observed average value in 
the local universe (Moustakas et al., 2006). The ratio may 
actually depend on the luminosity (e.g. i?-band in the 
rest- frame as in Moustakas et al. 2006), and adopting the 
fixed value may introduce some uncertainties. We need to 
examine the correlation between the Ha/ [OH] ratio and 
the luminosity (or metallicity) for star-forming galaxies 
at high redshift with extensive near-infrared spectroscopy 
as a future work. Figure 9 shows cumulative functions of 
SFR (left) and SSFR (right) of various samples. Blue lines 
and red lines correspond to field and cluster environments, 
respectively. We see an evolution of SFR function in the 
cluster environment in the sense that SFRs at highcr-z 
are systematically higher. However, SSFRs do not show a 
significant change over time. This indicates that the mass 
of galaxies which host star formation at z = 0.8 is smaller 
than that of z = 1.5. This is equivalent of the "down- 
sizing" effect (Cowie et al., 1996; Kodama et al., 2004; 
Thomas et al., 2005) where star- forming activity is shifted 
to low-mass galaxies with time. Also, there is a clear dif- 
ference between the two environments in the sense that 
SSFRs in the field environment are systematically higher 
than those in the clusters. 

In order to investigate the star formation activities at 
higher redshift, we make a comparison with the properties 
of galaxies in distant proto-clusters (Ha emitters at z = 2.2 
in PKS1138, Kurk et al. 2004; Ha emitters at z = 2.5 in 
4C23.56, Tanaka et al. submitted). Note however that 
the limiting SFRs for 4C23.56 only (- 25 Mgyr"!) is sig- 
nificantly larger than our limit for the other fields/clusters 
(12 M0yr^^), and the emitter samples in 4C23.56 do not 
include galaxies with lower SFRs. In Figure 10, we present 
the evolution of star formation activities both in the fields 
and in the clusters. Vertical axis shows the ratio of total 
star formation rate to total stellar mass, which is equiv- 
alent to the averaged specific star formation rate. There 
is a difference in the slope of evolution between the two 
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Fig. 9. The cumulative functions of SFR and SSFR. Red and 
blue lines show cluster and field samples, respectively: the 
Ho emitters a,t z = 0.8 in RXJ1716 (red dashed line; Koyama 
et al. 2010), [Oil] emitters at 2 = 1.2 in COSMOS field (blue 
solid line; Ideue et al. 2009), [Oil] emitters at z = 1.5 in 
XCS2215 (red solid line; Hayashi et al. 2010) and Ho emitters 
at 2 = 2.2 in GOODS-N field (blue dashed line; This work). 



environments. These two lines seem to cross at z ^ 2.5. 
This may indicate that, at high redshifts, we may reach to 
the epoch when star formation activities become higher in 
higher density regions than the general field, as opposed 
to the local relation. 

According to the Schmidt law (Kcnnicutt, 1998b), the 
SSFR can be related with the cold gas fraction, 

SSFR(x-P — . (5) 



where ^ = Mgas/ (Afgas + Mstar). In high density regions 
such as galaxy clusters, the gas fraction rapidly decreases 
from n ~ 0.4 at z = 2.5 to ^ 0.1 at z = 0.8. Here we de- 
termined the proportional constant by using the observed 
values of (SSFR) = 3 Gyr^^ and (/x) = 0.5 for the sample of 
Reddy et al. (2006). It implies that these galaxies actively 
form stars and consume a massive amount of gas on a 
shorter time-scale due to the environmental effect. In this 
argument, we have assumed a closed box model and ig- 
nored any merger effect which would affect the time-scale 
and complicate the argument. However, our main moti- 
vation here is to simply and roughly discuss the global 
evolution in SSFR and gas fraction. Note also that SSFR 
and gas fraction would not change if two identical galax- 
ies merge together. Therefore our discussion would not be 
significantly affected unless two galaxies at very different 
evolutionary stages merge together. 

The results presented in this paper on the time evo- 
lution and the environmental dependence of star form- 
ing activities are all intriguing. However, we are still 
very limited by the small area and/or the small number 
of observed fields and clusters. We are conducting the 
MAHALO-Subaru project (MApping HAlpha and Lines 
of Oxygen with Subaru; Kodama et al. in preparation) 
which aims to map out star forming activities over a 
wide redshift range (0.4 < z < 2.5) and in various envi- 
ronments by targeting many distant clusters and some 
general fields. By using a unique, large set of narrow- 
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Fig. 10. The evolution of star formation activities of the 
line emitting galaxies probed by narrow-band surveys. 
Filled circles and open circles indicate Ho emitters and 
[on] emitters, respectively. The 4C23.56 sample in- 
cludes only the galaxies with large SFRs (> 25M0yr-i). 



band filters available on the wide field cameras, Suprime- 
Cam and MOIRCS on Subaru Telescope, we can build a 
statistical, unbiased sample of line emitting galaxies (Ha 
and [OH]). This massive data set will enable us to sketch 
how the star formation activities propagate in the universe 
over the cosmic time and across environments, with much 
greater statistics and smaller uncertainties, hence provid- 
ing more comprehensive views of galaxy evolution in their 
most active phase of star formation and mass assembly. 

7. Summary 

We have conducted a narrow-band imaging and spectro- 
scopic surveys of Ha emitters at z = 2.2 in GOODS-North 
field, using MOIRCS on Subaru. Our survey has identified 
11 Ha emitters and one AGN over a 70 arcmin^ area. We 
have confirmed probable Ha emission lines for all of the 
seven targets by the spectroscopic follow-up observation, 
on top of the two already confirmed star forming galax- 
ies at z, 



spec 



2.2. Therefore, our technique of searching 
for Ha emitters at z = 2.2 based on the excess fluxes in 
the narrow-band (NB209) and photometric redshifts, is 
proven to be robust and efficient. The results and conclu- 
sions of this study arc summarized below: 

1. The Ha emitters have SFRs ranging from 
12 Mgyr^i to 60 Mgyr"!, with the mean SFR of 
(SFR) = 27.8]V[0yr-^ Note that we have corrected 
for dust extinction by assuming the typical value 
of A(Ha)=l. The averaged stellar mass of the Ha 
emitters is 4.0 x 10^" M©, and we find the correla- 
tion between the stellar mass and the specific star 
formation rate in the sense that more massive galax- 
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ies tend to have lower specific star formation rates. 

2. The Ha luminosity function is derived from our 
data by combining the data points reproduced from 
the previous works in the literature for the brighter 
magnitudes. The combined LF is represented by a 
Schechter function with logL = 42.82, log(/)= -2.78 
and a = —1.38. Our result shows a moderate steep- 
ness of the faint-end slope. By integrating the lu- 
minosity function thus derived, we find that the 
cosmic star formation rate at z = 2.2 is psFR = 
0.31 M0yr~^Mpc~^, which is consistent with other 
previous studies at z ^ 2. 

3. We have compared the properties of our emitters 
in the general field or low density environment, 
with those in the cluster environments to investi- 
gate the environmental dependence of galaxy evo- 
lution. There is a difference in the degree of time 
evolution of SSFR / SMgtar between the two en- 
vironments. This implies that the star formation 
activity is enhanced at z > 2 in high density regions 
as a consequence of "galaxy formation bias" in the 
early universe. 

We must warn however that there is a possibility that 
these results are still severely affected by the cosmic vari- 
ance, which may well be expected from the actual inho- 
mogeneous spatial distribution of our Ha emitters. It is 
necessary to extend the survey area to cover a represen- 
tative volume of the universe and average over the cosmic 
variance in order to obtain more robust conclusions. 

Acknowledgments 

This paper is based on data collected at Subaru 
Telescope, which is operated by the National 
Astronomical Observatory of Japan. We thank the 
Subaru telescope staff for their help in the observation. 
We are also very grateful to the MODS team for allowing 
us to use MODS data and catalog and Yuko Ideue for 
providing information about [OH] emitters at z ~ 1.2 
in COSMOS field. We thank the anonymous referee 
who gave us many useful comments, which improved 
the paper. T.K. acknowledges the financial support in 
part by a Grant-in-Aid for the Scientific Research (Nos. 
18684004 and 21340045) by the Japanese Ministry of 
Education, Culture, Sports, Science and Technology. 
Y.K. acknowledge the support from the Japan Society for 
the Promotion of Science (JSPS) through JSPS research 
fellowships for young scientists. 

References 

Alexander D. M., et al. 2003, AJ, 126, 539 
Arnaboldi M., et al. 2003, AJ, 125, 514 
Atek H., et al. 2010, ApJ, 723, 104 
Bertin E., Arnouts S., 1996, A&AS, 117, 393 
Bruzual C, Chariot S., 2003, MNRAS, 344, 1000 
Bunker A. J., Warren S. J., Hewett P. C, Clements D. L., 
1995, MNRAS, 273, 513 



Calzetti D., Armus L., Bohlin R. C, Kinney A. L., 
Koornneef J., Storchi-Bcrgmann T., 2000, ApJ, 533, 
682 

Chabrier C, 2003, PASP, 115, 763 

Chapman S. C, Blain A. W., Small I., Ivison R. J., 2005, 
ApJ, 622, 772 

Cowie L. L., Songaila A., Hu E. M., Cohen J. C, 1996, 
AJ, 112, 839 

Daddi E., Cimatti A., Renzini A., Fontana A., Mignoli 
M., Pozzetti L., Tozzi P., Zamorani C, 2004, ApJ, 617, 
746 

Dale D. A., et al. AJ, 135, 1412 

Doherty M., Bunker A., Sharp R., Dalton C, Parry I., 

Lewis I., 2006, MNRAS, 370, 331 
Doherty M., et al. 2004, MNRAS, 354, L7 
Elbaz D., et al. 2007, A&A, 468, 33 

Erb D. K., Steidel C. C, Shapley A. E., Pettini M., Reddy 

N. A., Adelberger K. L., 2006a, ApJ, 647, 128 
Erb D. K., Steidel C. C, Shapley A. E., Pettini M., Reddy 

N. A., Adelberger K. L., 2006b, ApJ, 646, 107 
Fujita S. S., et al. 2003, ApJL, 586, L115 
Gallego J., Zamorano J., Aragon-Salamanca A., Rego M., 

1995, ApJL, 455, Ll-f 
Geach J. E., Small I., Best P. N., Kurk J., Casah M., 

Ivison R. J., Coppin K., 2008, MNRAS, 388, 1473 
Hayashi M., Kodama T., Koyama Y., Tanaka I., 

Shimasaku K., Okamura S., 2010, MNRAS, 402, 1980 
Hayes M., Schaerer D., Ostlin G., 2010, A&A, 509, L5-h 
Hippelein H., et al. 2003, A&A, 402, 65 
Hopkins A. M., 2004, ApJ, 615, 209 
Hopkins A. M., Beacom J. F., 2006, ApJ, 651, 142 
Ideue Y., et al. 2009, ApJ, 700, 971 
lye M., et al. 2004, PASJ, 56, 381 
Kajisawa M.,et al. 2009, ApJ, 702, 1393 
Kennicutt Jr. R. C, 1998a, ARA&A, 36, 189 
Kennicutt Jr. R. C, 1998b, ApJ, 498, 541 
Kennicutt Jr. R. C, Tamblyn P., Congdon C. E., 1994, 

ApJ, 435, 22 

Kinney A. L., Calzetti D., Bohlin R. C, McQuadc K., 
Storchi-Bergmann T., Schmitt H. R., 1996, ApJ, 467, 
38 

Kodama T., Arimoto N., Barger A. J., Arag'on- 

Salamanca A., 1998, A&A, 334, 99 
Kodama T., BeU E. F., Bower R. G., 1999, MNRAS, 302, 

152 

Kodama T., et al. 2004, MNRAS, 350, 1005 

Koyama Y., Kodama T., Shimasaku K., Hayashi M., 

Okamura S., Tanaka I., Tokoku C, 2010, MNRAS, 403, 

1611 

Koyama Y., et al. 2008, MNRAS, 391, 1758 
Ly C, et al. 2007, ApJ, 657, 738 

Moorwood A. F. M., van der Werf P. P., Cuby J. G., Oliva 

E., 2000, A&A, 362, 9 
Morioka T., Nakajima A., Taniguchi Y., Shioya Y., 

Murayama T., Sasaki S. S., 2008, PASJ, 60, 1219 
Moustakas J., Kennicutt Jr. R. C, Tremonti C. A., 2006, 

ApJ, 642, 775 



Cosmic Star Formation Activity at 



Nakagawa T., 2008, in Society of Photo-Optical 
Instrumentation Engineers (SPIE) Conference Series 
Vol. 7010 of Society of Photo-Optical Instrumentation 
Engineers (SPIE) Conference Series, SPICA mission for 
mid- and far-infrared astronomy 

Ouchi M., et al. 2004, ApJ, 611, 660 

Perez-Gonzalez P. C, Zamorano J., Gallego J., Aragon- 
Salamanca A., Gil de Paz A., 2003, ApJ, 591, 827 

Pickles A. J., 1998, PASP, 110, 863 

Reddy N. A., Steidel C. C, Erb D. K., Shapley A. E., 

Pettini M., 2006, ApJ, 653, 1004 
Reddy N. A., Steidel C. C, Fadda D., Yan L., Pettini M., 

Shapley A. E., Erb D. K., Adelberger K. L., 2006, ApJ, 

644, 792 

Rousselot P., Lidman C, Cuby J., Moreels C, Monnet 

G., 2000, A&A, 354, 1134 
Salpeter E. E., 1955, ApJ, 121, 161 
Schcchtcr P., 1976, ApJ, 203, 297 
Shioya Y., et al. 2008, ApJS, 175, 128 
Skrutskie M. F., et al. 2006, AJ, 131, 1163 
Sobral D., et al. 2009, MNRAS, 398, 75 
Steidel C. C, Shapley A. E., Pettini M., Adelberger K. L., 

Erb D. K., Reddy N. A., Hunt M. P., 2004, ApJ, 604, 

534 

Suzuki R., et al. 2008, PASJ, 60, 1347 

Tanaka M., De Breuck C, Venemans B., Kurk J., 2010, 

A&A, 518, A18+ 
Teplitz H. I., Malkan M., McLean I. S., 1998, ApJ, 506, 

519 

Thomas D., Maraston C, Bender R., Mendes de Oliveira 

C, 2005, ApJ, 621, 673 
Thompson D., Djorgovski S., Beckwith S. V. W., 1994, 

AJ, 107, 1 

Tokoku C.,et al. 2006, in Society of Photo-Optical 
Instrumentation Engineers (SPIE) Conference Series 
Vol. 6269 of Society of Photo-Optical Instrumentation 
Engineers (SPIE) Conference Series, Infrared multi- 
object spectrograph of MOIRCS 

Tran K., et al. 2010, ApJL, 719, L126 

Villar v., Gallego J., Perez-Gonzalez P. G., Pascual S., 
Noeske K., Koo D. C, Barro G., Zamorano J., 2008, 
ApJ, 677, 169 

Yagi M., Kashikawa N., Sekiguchi M., Doi M., Yasuda N., 
Shimasaku K., Okamura S., 2002, AJ, 123, 66 

Yan L., McCarthy P. J., Freudling W., Teplitz H. I., 
Malumuth E. M., Weymann R. J., Malkan M. A., 1999, 
ApJL, 519, L47 



